Rheology of colloidal microphases in a model with competing interactions 
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We study the rheological properties of colloidal microphases in two dimensions simulating a model 
of colloidal particles with competing interactions. Due to the competition between short-range 
attraction and long-range repulsion, as a function of the density the model exhibits a variety of 
microphases such as clusters, stripes or crystals with bubbles. We prepare the system in a confined 
microphase employing Monte-Carlo simulations and then quench the system at T = 0. The resulting 
configurations are then sheared by applying a drag force profile. We integrate numerically the 
equation of motion for the particles and analyze the dynamics as a function of the density and 
the applied strain rate. We measure the stress-strain curves and characterize the yielding of the 
colloidal microphases. The results depend on the type of microphase: (i) clusters are easily sheared 
along layers and the relative motion is assisted by rotations, (ii) Stripes shear easily when they are 
parallel to the flow and tend to jam when are perpendicular to it. Under a sufficiently strong shear 
rate perpendicular stripes orient in the flow direction, (iii) Crystals with bubbles yield by fracturing 
along the bubbles and eventually forming stripes. We discuss the role of dislocations, emitted by 
the bubbles, in the yielding process. 



I. INTRODUCTION 

Examples of spontaneous pattern formation or micro- 
separation, in condition close to the thermodynamic equi- 
librium, can be found in many different experimental sys- 
tems, both in two dimensions [l], H, H, El and in three 
dimensions [|| 0, H, • The interest in these novel phases 
of matter is very high because the skill in controlling the 
architecture of particle aggregates, as well as the super- 
structures the latter arrange on, is nowadays a strategic 
tool in nanotcchnology. in order to engineer new materi- 
als with specific electric, magnetic, optical and rheologi- 
cal properties. 

The topic of pattern formation is important not only 
from a technological point of view, but also from a fun- 
damental one, influencing many different branches of sci- 
ence: from vitrification and gelationpjj [ll], Ell [U| to 
colloidal systems [Til, flfil llH Uj\ , from biological mem- 
branes [3 [l9[ to network formation of tissue cells [2(| , 
to physisorbed layers on solid surfaces [2l|, • Even if 
all of these systems are rather different at the molecular 
level, the patterns formed by the particles can be sim- 
ilar: circular clusters, stripes and rings of particles can 
be easily observed in most of the bidimensional systems, 
while lamellae, spherical clusters and cylinders occur in 
three dimensions. Moreover the particle aggregates can 
be often arranged to form disordered liquid-like configu- 
rations, but also crystalline-like super-structures such as 
grids of stripes and lamellae, hexagonal lattice of clusters 
and bubbles. 

In many cases a simple mechanism to explain the 
micro-separation is the competition between a short- 
range attraction and a long-range repulsion. For col- 
loidal systems, the short-range attraction can be due 



to depletion forces and Van der Waals forces, while the 
longer-range repulsion stems from dipolar forces as in 
Langmuir monolayers and ferrofluids, or from partially 
screened electrostatic forces. For sufficiently strong long- 
range repulsion, the standard vapor to liquid transi- 
tion is inhibited in favor of the microseparation. In 
particular, the short-range attraction favors the con- 
densation of particles, while the growth of the clusters 
is limited by the long-range repulsion. Recently how- 
ever, models based on purely repulsive interactions have 
been discussed which also support the micro-separations 

muni. 

Although great attention has been given to the study 
of the micro-separated systems via theoretical llCt l26l . 
l27l [IH and numerical approaches [29|, H3, HH, |32| , l33l |. 
a general comprehension of the behavior of such phases 
is still lacking. Recently much attention has been also 
given to the role of confinement on pattern formation: 
simulations and experiments suggest, in fact, that con- 
finement is a powerful tool to modify and even induce 
new pattern morphologies, which in the bulk might not 
appear at all. Most of these studies concern 3D block 
copolymer systems between smooth or patterned sub- 
strates [33 . [3fil l36l l37ll and recently 2D colloidal systems 
under confinement |38| |. 

Many common complex fluids such as pastes, emul- 
sions, foams, colloidal systems and granular materials are 
studied also from the point view of the mechanical and 
rheological properties, which are relevant for industry. 
These systems, in fact, might exhibit both solid-like and 
liquid-like responses depending on the length scale, the 
time scale and on the applied shear stress. Moreover the 
shear has proved to be a mechanism capable of inducing 
and inhibiting pattern formation as well. Many experi- 
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ments involve aqueous foams under shear, which act as 
elastic solids for small deformations, while flow for large 
applied shear [H, H(|. Moreover cluster phases in col- 
loidal systems, which become dilatant under flow, have 
been recently reported in (4ll . |42| . Finally in Ref. (43| 
the role of geometrical confinement in the shearing of 
colloidal systems is considered. 

Simulations in 3D have pointed out that shear can fa- 
vor or inhibit the cluster phase with respect to the homo- 
geneous one, via shearing [l4|- In particular in Ref. [44|, 
the authors have analyzed the behavior of mixtures of 
particles both in the brownian regime, where the shear 
forces are comparable to the brownian ones, and in the 
Stokes regime where the shear forces greatly exceed the 
brownian forces. In both cases, systems which form gels 
at rest, support the formation of much more compact 
and somehow complex structures under shear. The role 
of external stresses, close to the jamming transition and 
during the aggregation processes in attractive colloidal 
systems, is presented in |4j|, where a general phase di- 
agram is proposed as a function of density temperature 
and stress. 

Concerning rheological properties of 2D systems, most 
of the experiments focus on the homogeneous phases of 
Langmuir monolayers 0, [13]. Recently attention has 
been paid also to the inhomogeneous phases as in Ref. 

, where a single layer of bubbles floating on a liquid 
surface is studied. Very interestingly such a system shows 
the coexistence between a flowing and a jammed phase of 
bubbles. Simulations instead address mainly the proper- 
ties of foams as in [49L l50l . l5l| , where flow is strongly in- 
termittent, once the shear stress exceeds the yield stress. 
The avalanche-like rearrangement of the bubbles are vis- 
ible from the stress-strain curve, where the fluctuations 
are strongly enhanced. In other words the stress is at first 
accumulated by macroscopic deformation of the struc- 
ture and then abruptly released, followed by topological 
changes. Stain bursts are typical in the plastic deforma- 
tion of amorphous solids, as shown in numerical simula- 
tions [HI, [H, [H| , but also of microscale crystals [HI, [56| 
where they are due to the collective dynamics of disloca- 
tions [13, Eli ■ 

In this paper we present results of numerical simula- 
tions on the rhcology of two dimensional colloidal mi- 
crophases in confined geometries. We consider an ensem- 
ble of colloidal particles interacting by a potential with 
short-range attraction and long-range repulsion. As a 
function of the density and depending on the constraints 
imposed by the confinement, the colloidal system displays 
a variety of equilibrium microphases, such as clusters, 
stripes or bubbles [38|. Here, we analyze the behavior of 
the system when a shear profile is imposed on the fluid. 
In particular, we focus on the stability of the pattern 
morphology under shear, especially in the intermediate 
high-density region which has not been studied so far. 
The rheological behavior of the system is dependent on 



the density and on the corresponding microphase. At 
low densities, the colloidal particles separate into clus- 
ters that are easily sheared. In this case we observe 
that shearing is assisted by cluster rotations, a process 
that is reminiscent of the shear of granular media [59| . 
At higher densities the colloidal particles form stripes. 
When the stripes are oriented parallel to the shear di- 
rection, the system deforms easily, following the imposed 
velocity profile. On the other hand, when the stripes are 
oriented perpendicular to the shear direction the system 
tends to jam. Only under the application of sufficiently 
high shear rate the system is able to flow. This is done by 
a reorientation of the stripes along the shear direction. 
Similarly, bubble phases are jammed at low shear rates 
and flow at large shear rate. In order to flow, the bub- 
bles tend to coalesce into stripes which are then easily 
sheared. 

The aforementioned rheological behavior is reflected in 
the stress strain characteristics of the system. In the case 
of perpendicular stripes and bubbles, the shear stress dis- 
plays a yield point, corresponding to the break-up of the 
jammed configurations, followed by a shear weakening 
regime. For large shear strain, the stress reaches a con- 
stant value corresponding to the flow of parallel stripes. 
Shearing a system which is prone to microseparate is 
similar to the application an external force field which 
favors specific pattern morphologies or specific orienta- 
tions. The possibility to control, modify and stabilize 
a specific pattern is particularly appealing for the tech- 
nology of new materials. From this point of view, the 
application of shear stress could represent a promising 
strategy to control the structure of microphases. 

This paper is organized as follows. First we describe 
the microscopic model of the system and its equilibrium 
properties, in the bulk and under confinement, in sec- 
tion III At then the rheological quantities as well as the 
simulations details are introduced in section III Bl Then 
in sections IIIII HVl and |Vl respectively, we describe the 
rheological properties for the cluster, the stripe and the 
bubble phases. In particular, for each case we discuss 
the shear rate profile as a function of time, the flow curve 
(that is the stress-strain curve) and the velocity profile of 
the particle aggregate as a function of the distance from 
the confining walls. A comparison between the behav- 
ior of the fluid under shear, adopting different pattern 
morphology, is drawn in section I VII 



II. THE MODEL 

A. Equilibrium 

Within the framework of the competing interactions, 
we consider a bidimensional fluid, in which the effective 
potential for a pair of colloidal particles is: 



3 



U (r) = \ A W ex P(-£) + %■ ex P("ir) if r ^ (i) 

I otherwise 



r being the intcrparticle distance and a the particle 
diameter. Such a potential is characterized by a soft-core 
repulsive term at very short distances, followed by an 
attractive well plus a repulsive hump at larger distances 
(see Fig. [1] upper panel). The potential parameters 
are: n = 12, A = 0.018 for the short-range repulsion, 
R a = 1 <T,e a = 1 for the short-range attraction and 
R r = 2 a, e r = 1 for the longer-range repulsion. The 
short-range repulsion takes into account the impene- 
trability of the particles, while the remaining part of 
the potential mimicks the colloid interaction mediated 
by the surrounding fluid. The potential is truncated 
at R cu t = 10 a. With the present choice of the 
parameters, the minimum of the potential is located at 
r ~ 1.2cr, while the maximum value of the longer-range 
repulsive hump is at r ~ 4er . 
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The interaction potential of equation ([T]) is similar to 
those studied in [TO, [H, H3 f° r three dimensional sys- 
tems. In case of two-dimensional fluids our model is 
the same discussed in [3l|, [38|, but for the short-range 
repulsion which is treated as a hard-core in the latter, 
while it is implemented as a soft-core potential in the 
present work. Such a difference is not particularly im- 
portant from the point of view of the general features 
of the phase diagram in bulk, which is schematically de- 
picted in Fig. [T] (bottom panel). At sufficiently low tem- 
perature, the present model supports the passage from 
droplets to stripes to bubbles as the density increases. 



homogeneous fluid phase 



The main difference with the bulk model studied in 31 1 
is merely an overall shift of the phase diagram towards 
lower densities. 

In our model, the fluid is confined along the x direction 
by smooth parallel walls, the separation of which is L x . 
The wall-particle interaction is treated as a rather steep 
potential of the form 



U wp (x) = C exp(-x/£), 



(2) 



where x is the distance from the wall. The parameters 
we use in simulations are £ = 0.1c and C = Uq, Uq being 
the depth of the attractive well of the particle-particle 
potential. As the range of the wall-particle interaction is 
very short and the repulsion very strong, such a poten- 
tial mimicks very well the effects of neutral hard walls, 
as we have also verified performing some Monte Carlo 
simulations with a hard-wall potential too. 

Under lateral confinement, if the hard disk potential 
is used for the short-range repulsion as in [111, a transi- 
tion from stripe to clusters occurs as the temperature is 
lowered. Such a transition disappears in our case, where 
the soft-core potential repulsion is used. The results dis- 




FIG. 1: Upper panel: particle-particle effective potential corre- 
sponding to Eq. {T] in units of the potential well depth Uo = 
\U pp (r = cr) | ) . Bottom panel: schematic phase diagram, in the 
density, temperature plane (p, T) for a bulk bidimcnsional fluid 
subject to competing interactions. 



cussed in (38| can be obtained if the exponent n of Eq. ([T]) 
is increased up to 36. But as the aim of this work is study- 
ing the effect of shear on different pattern morphologies, 
disregarding at first the thermal noise, the change of mor- 
phology with temperature is not of interest at this stage. 
Moreover, with the present choice of the potential param- 
eters, we are able to study the stripe and bubble phase 
at relatively lower densities, meaning that we carry out 
simulations with smaller systems. 

In this work we discuss three densities 
pa 2 = 0.15, 0.3, 0.5, corresponding to the clus- 
ter, stripe and bubble regime. The clusters and the 
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bubbles are arranged onto a triangular super-lattice, so 
that the particle configurations are specular. One of 
the goals of our work, in fact, is to study how such a 
hypothetical symmetry between filled and empty region 
in the particle configuration affect the rheology of the 
fluid. The stripes, instead, form a grid of parallel layers. 
In each microphase the pattern period is P ~ 11 a. The 
crystalline order of the super-lattice is signaled by the 
occurrence of Bragg peaks in the static structure factor 
at short wave vectors whose length is roughly connected 
to the period through the relation P Kt 2ir/K p . Random 
phase approximation for this kind of potential has 
provided predictions in good agreement with simulations 
[3lj . The maximum size for the particles aggregates, 
instead, is dominated by the width of the attractive well 
of the particle-particle potential. So that, in the present 
case, we see that cluster diameter is ~ 4 er, similarly to 
the bubble diameter and to the stripe width. 

The addition of a geometrical constraint, such as the 
presence of neutral walls along one of the directions, is 
a source of further frustration for the system, mainly if 
the wall separation L x is not commensurate to the in- 
trinsic bulk periodicity. For our model, the most pecu- 
liar case is represented by the striped phase, where the 
alignment of the stripes depends strongly on the value 
of L x . In particular, if L x — mP + A s (m is an inte- 
ger and A s the stripe width), the stripe orientation is 
parallel to the walls, otherwise we observe the forma- 
tion of a mixed phase with two parallel stripes next to 
the walls and stripes perpendicular to the walls in the 
center of the slit. In figure [2] typical configurations for 
different patterns between walls are plotted. Such config- 
urations are generated through Monte Carlo simulations 
in the canonical ensemble (NVT, N number of particles, 
V volume, T temperature). In our simulations N ranges 
from 1000 to 2000 particles, and the temperature is equal 
to T = 0AU o . 



B. Rheology 

To simulate the rheology of the colloidal system, we 
consider a simple overdamped dynamics for the parti- 
cles under the action of an external force field due to the 
shear. A moving colloidal particle in a fluid is subject 
to a drag force. In particular, for a spherical particle 
of radius a moving at small velocity v, the drag force is 
given by Fd = —Vv, with T = Qirrja where r\ is the fluid 
viscosity. In this paper, we apply a weak linear shear 
profile to the carrier fluid V(x) — Voxy/L x and, corre- 
spondingly, the colloid particle is subject to a drag force 
profile Fd{x) = Foxy/L x , where Fq — TVq. Neglecting 
inertial and thermal effects, the equation of motion for 
the particles is given by: 
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FIG. 2: Typical configurations obtained via Monte Carlo simula- 
tions for different densities and different values of the wall separa- 
tion L x and T = OAUo. We notice that when the wall separation 
L x is not commensurate to the stripe period, the perpendicular 
alignment is favored upon the parallel one. 



where the interparticle forces are computed from the po- 
tential in Eq. ([T]) 



fpp(r) 



dU p 
df 



(4) 



Eq. ([3]) is integrated numerically by means of an adaptive 
step-size fourth-order Runge-Kutta algorithm. In the 
present simulations, we fix the timescales setting T = 1. 
Equilibrium configurations obtained from Monte-Carlo 
simulations at T = OAUo are first relaxed at T = in 
absence of external forces, integrating Eq. ([3]) for Fq = 0. 
Then we apply an external shear profile to the resulting 
configurations and study the system as a function of F 
and p. 

It will be instructive to compare the internal shear 
strain rate with the externally imposed strain rate j ex t = 
\^Sx = Vo/(2L x ). The internal strain rate is in general 
not uniform and it is necessary to consider its average 
value defined as 



<7> 



[<v+>-<v y >] , 



(5) 



n), 



(3) 



where < > is the y component of the particle velocity 
restricted to the x ^ domain. To obtain a clearer 
picture of the internal strain rate, in order to quantify the 
deviations from a laminar flow, we measure the particle 
velocity profile (v y (x)) dividing the system into n s strips 
oriented along y and averaging the y component of the 
velocities of the particles in each strip. In particular, 
we use n s = 7 for L x = 76 and n s = 5 for L x — 48, 53. 



Finally, wc measure the stress-strain curves of the system 
computing the internal shear stress as 

Vxy = -Y~r J2( Xi ~ x j)fppH^ ~ ( 6 ) 

Ll X Lly 

» I] 

III. CLUSTER PHASE 

The rheological properties of the cluster phase are 
ruled by the interactions between individual clusters. At 
low temperatures individual particles inside the clusters 
do not move, thus the clusters can be considered as rigid 
objects which repel each other. In Fig. [3] (upper panel) 
we report the average strain rate 7 of the system normal- 
ized by the external strain rate. The average strain rate 
shows characteristic oscillations, which are a peculiarity 
of the pattern geometry: at low temperatures the clus- 
ters arrange into a regular super-lattice and to flow the 
system must overcome a Peierls-Nabarro periodic poten- 
tial. In other words, the oscillations are connected to the 
sliding past of clusters moving along adjacent columns. 
The period of the oscillations diminishes as the external 
applied force Fq increases, as the clusters move faster. 
Oscillations are also present in the stress-strain curves as 
shown in Fig. [3] (bottom panel). For the timescales and 
shear rates that we have studied, we see a slowing down 
of the dynamics at low Fq but the system does not jam. 
We would expect that at very low Fq the system could be 
blocked by the weak intra-cluster potential, but we could 
not simulate efficiently a similar regime as the relaxation 
time becomes too long. 

In Fig. 0] (upper panel) we show the velocity profile of 
the clusters. The external applied velocity profile is also 
shown as a dashed line. At low Fq the velocity profile 
shows that the colloid particles are slower with respect to 
the surrounding medium. Increasing the external force, 
instead, the colloid velocity profile tends to be equal to 
the external one (dashed line). In particular, as Fq in- 
creases, the velocity profile tends to be linear at first in 
the middle of the box and then nearby the walls. Hence 
the flow is at first newtonian in the middle of the box, 
while slippage can be present close to the walls. Such 
a behavior can be considered as an indication of shear 
banding, where different spatial regions of the fluids are 
characterized by different flowing properties such as dif- 
ferent values of the viscosity. 

We see that the cluster lattice is distorted by the shear 
deformation and clusters in different columns have to 
slide past each other. This process is occasionally helped 
by rotations around the clusters center of mass. Due to 
the shear profile individual clusters are subject to a net 
torque that can induce rotations, creating a sort of roller- 
bearing effect (See EPAPS Document No. XXXX.l for 
an animation of the shear of the cluster phase). Similar 
rotations are commonly observed in the shear of gran- 
ular media and is very important in earthquakes when 
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FIG. 3: (Color online) Upper panel: the average strain rate, 
divided by the external strain rate ■jext = Va/ (2L X ), as a function 
of time for different values of Fq for a system with N = 1000, 
p = 0.15 and L x = 76. Bottom panel: the stress-strain curves for 
the same cases shown in the upper panel. 



the faults are filled by a granular gouge (59j . Similarities 
between colloidal systems and granular media have been 
pointed out on the basis of rheological measurements in 
Ref. [4l[. Here, we measure the average torque of the 
clusters as a function of the horizontal coordinate of their 
center of mass. As can be seen in Fig. @] (bottom panel) 
the torque displays a characteristic pattern depending 
on the cluster position. The clusters nearby the walls 
rotate less than the rest because the interaction with the 
surrounding clusters is limited on one side; the central 
clusters, instead, are those for which the applied velocity 
is less important as most of the particles experiment an 
almost zero drag velocity, so that their rotation is essen- 
tially due to the repulsion with the clusters flowing in 
the nearby columns. Finally the clusters whose position 
is intermediate between the walls and the center of the 
simulation box, experiment both the effects due to the 
applied drag force and to the cluster-cluster repulsion, so 
that the net torque is the highest. 
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Fig. [5l upper panel). This behavior is reflected by the 
stress strain curve (Fig. [51 bottom panel), that changes 
character for Fq > 0.03 where the yield stress is followed 
by a strain-rate weakening part while, for Fq < 0.03 
the stress reaches a plateau. The peak stress for low 
Fq scales in a viscous manner as dp ~ Cj and changes 
character for Fq > 0.03 scaling as a p ~ ay + C'j, with 
C < C (inset of Fig. [S]). We also note strong fluctua- 
tions in the high strain regime, a behavior which is of- 
ten observed in systems with avalanche-like flows as in 

0,m[H,[5l[5l[5l,[63. 

The jamming of the system is clearly visible in the 

velocity profiles, where we see that the central part of 
the system has zero velocity and only the two boundary 
parallel stripes arc flowing (Fig. [5J upper panel) . Typical 
snapshots of the system, when it is jammed and when it 
is flowing are plotted in the bottom panels of Fig. [6] too 
(See EPAPS Document No. XXXX.2 for an animation of 
the shear of the perpendicular stripe phase showing the 
reorientation of the stripes at the yield point). 



The case of parallel stripes is particularly simple to 
shear. The velocity profile rescaled to the applied force 
is always linear for each value of Fq we have used. The 
most interesting aspect is perhaps the behavior of the 
stress strain-curves which shows some fluctuations, even 
if much more smaller than those observed with perpen- 
dicular stripes, and probably due to the irregularity in 
the stripe surfaces (see Fig. [7]). The mean stress is a lin- 
ear function of the strain rate, indicative of a simple vis- 
cous behavior. An extrapolation at low shear rate would 
indicate a very low yield stress of the order of ay — 10~ 5 . 



FIG. 4: (Color online) Upper panel: the velocity profiles for 
different values of Fq, for a system with N = 1000, per 2 = 0.15 
and L x = 76er. Notice the deviation from the externally imposed 
linear profile (dashed line) occurring at low shear rates. Data are 
averaged over time and over 10 realizations. Bottom panel: the 
torque M profiles acting on the clusters. In the inset we show a 
typical configuration of the system. 



IV. STRIPE PHASES 



When we consider the rheology of the stripe phase, we 
have to pay attention to the orientation of the stripes 
with respect to the shear direction. By slightly tuning 
L x we can obtain configurations in which the stripes are 
all parallel to the walls or mixed configurations made up 
of two stripes parallel to the walls, while in the middle 
they arrange perpendicularly to the walls. 
The most interesting case is that corresponding to the 
perpendicular stripes. In fact, at low forces the system 
jams, failing to reorient in the flow direction. Hence, for 
Fq < 0.03, the average strain rate is lower than the exter- 
nal rate, while for larger Fq their ratio tends to one (see 





x/L 




FIG. 5: (Color online) Upper panel: the average strain rate as a 
function of time for different values of Fo for a system with N = 
1000, p = 0.3 and L x = 53 plotted in log-linear scale. In the initial 
condition the stripes are mostly perpendicular to the flow direction. 
Bottom panel: the stress-strain curves. In the inset we report the 
maximum stress a v as a function of the strain rate 7. Notice the 
change of behavior as the systems yields at higher rates. 



FIG. 6: (Color online) Upper panel: the velocity profiles as a func- 
tion of time for different values of i*o, for a system with TV = 1000, 
p = 0.3 and L x = 53. At low shear rate only the outer (paral- 
lel) stripes moves while the perpendicolar stripes are jammed. At 
higher shear rate the perpendicolar stripes break and flow. The 
dashed line represents the imposed profile. Bottom panels: parti- 
cle configurations under different conditions. Initial configuration 
(left); yield point, when the central stripes break at F = 0.1 and 
7 = (middle); high strain Fg = and 7=1 (right). 
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FIG. 7: (Color online) The stress-strain curves for different values 
of F for a system with N = 1000, p = 0.3 and L x = 48. The 
initial state is characterized by stripes which arc all parallel to the 
walls and to the flow direction. 
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V. BUBBLE PHASES 

Increasing the density, the system develops a bubble 
phase in which the bubbles are arranged on a triangu- 
lar superlattice. From a qualitatively point of view, the 
bubbles under shear become more and more elongated 
along the direction of the applied velocity field, until they 
join to form stripe-like configurations (bottom panel of 
FiglH]). The rheology of the bubble phase has some sim- 
ilarities with the one of the perpendicular stripe phase, 
since both phases jam at low shear rates. The main dif- 
ference with respect to the stripe phase is that for small 
external applied strain, the bubble phase is completely 
jammed. This is evident from the slowing down of the 
average strain rate for F < 0.005 (Fig. [8] upper panel) 
and also from the velocity profile which is completely 
flat (Fig. upper panel). In the case of perpendicular 
stripes, instead, the jamming involves only the central 
region, while the stripes close to the walls flow as well. 

When the applied force to the bubble phase is greater 
than the peak stress, the system starts to flow and once 
again parallel stripes arc formed (Fig. [5] bottom panel), 
since this state is the easiest to shear (See EPAPS Docu- 
ment No. XXXX.3 for an animation of the shear induced 
stripe formation in the bubble phase). The peak stress is 
an increasing function of the strain rate as depicted in the 
inset of the bottom part of Fig. [8] Extrapolations at low 
shear rate suggest an yield stress of the order ay = 0.02. 
The shear weakening observed in the stress-strain curve 
is reminiscent of the stripe phase, but the yield stress 
cry measured in the bubble phase is about an order of 
magnitude greater than that observed in the stripe case. 
This is easily understandable as, in the bubble phase, the 
particles are strongly packed together, so that the motion 
of the individual particles is hindered by the surronding 
ones. In other words, collective motion of particles is nec- 
essary to ensure a structural change of the system. We 
notice as well that the stress-strain curve is characterized 
by large fluctuations, suggesting again intermittent and 
abrupt rearrangements as commonly observed in plastic- 
ity [51 HE [H m an d in foam rheology [IE [IE ED] ■ 

In order to better understand the yielding of the bub- 
ble phase, we perform a voronoi triangulation of the par- 
ticle system and follow the evolution of the topological 
defects. Due to the T = condition, the system is locally 
ordered with a few topological defects around the bub- 
bles. As the system is sheared, we observe the nucleation 
and propagation of dislocations, corresponding to pairs of 
five-fold and seven-fold coordinated atoms (See EPAPS 
Document No. XXXX.4 for an animation of the voronoi 
triangulations of the configurations in the sheared bubble 
phase). I n particular, dislocations are sometimes created 
at the surface of a bubble and then propagate towards 
the nearest bubble (see Fig. [TUJ). Hence, we can see the 
shear induced transformation from bubbles to stripes as a 
ductile fracture process: the system first deforms plasti- 
cally in the vicinity of the bubbles and eventually a crack 
propagates, leading to the coalescence of the bubbles into 



stripes. 

It is also interesting to compare the rheology of the 
bubble phase with the one of the cluster phase. At equi- 
librium, the cluster phase and the bubble are character- 
ized by configurations which are specular. Under shear, 
however, the clusters are easily sheared also for very low 
values of Fq, while a yield stress must be overcome in the 
bubble phase. Moreover, while the clusters rotate and 
mantain their shape until Fo is rather large, the bub- 
bles deform almost immediately. Notwithstanding these 
large differences in the structure of the systems, the ve- 
locity profiles appear to be very similar for Fq > 0.005. 
In particular the velocity profiles obtained in the bubble 
phase are similar to those obtained in the cluster phase 
but with smaller external forces. Moreover, the veloc- 
ity profiles tend to become linear at first in the middle 
of the box simulation and eventually nearby the walls, 
indicating once again the presence of shear banding. 




FIG. 8: (Color online)Upper panel: the average strain rate as 
a function of time for different values of Fo for a system with 
N = 2000, per 2 = 0.5 and L x = 76 a plotted in log-linear 
scale. Notice that at large Fo the strain rate reaches a steady 
value, while for Fo = 0.0001 the strain rate decreases towards zero, 
indicating a jammed phase. Bottom panel: the stress-strain curves 
for different values of Fo- In the inset wc show the peak stress as 
a function of the applied strain rate 7. 




FIG. 9: (Color online) Upper panel: the velocity profiles for dif- 
ferent values of the drag force Fo for a system with N = 2000, 
per 2 = 0.5 and L x = 76a. Data are averaged over time and over 
10 realizations. Notice the deviation from the externally imposed 
linear profile occurring at lower shear rates. For Fo = 0.0001, the 
system is jammed and the profile is flat. Bottom panels: subse- 
quent snapshots for the system at Fo = 0.1. The snapshots are 
taken at the beginning of the simulation (left), at the yield strain 
(middle) and at large strain (right). Notice the crossover from 
bubble to stripes at large strains. 



FIG. 10: (Color online) Voronoi triangulations ol the config- 
uration ol the system configurations at two subsequent time 
steps for a system with N = 2000, p = 0.5, L x — 76 and 
Fo — 0.1. The snapshot are taken around the yield strain 
7 ~ 0.11. Pairs of seven-fold and five-fold particles form dis- 
locations and are colored respectively in green and red. For 
clarity, we do not show the topological defects arising at the 
sample boundary and around the bubbles. Notice the motion 
of a dislocation between neighboring bubbles (indicated by an 
arrow) . 



VI. DISCUSSION AND CONCLUSIONS 

In simple liquids the flow properties are entirely de- 
termined by the viscosity which does not depend on the 
shear rate. In the laminar regime, for simple liquids a 
linear relation between a and 7 holds. For the system 
we study, such behavior is not present anymore, because 
particle aggregates lead to a coupling between the mi- 
crostructure of the system and its flow. Sometimes, the 
coupling between microstructure and flow is so strong 
that we obtain pattern morphologies which are not pos- 
sible at equilibrium for a fixed density or for a fixed wall 
separation. 

To summarize the behavior of the system at different den- 
sities, wc compare the strain rate obtained for different 
pattern morphology. In particular, the strain rate is av- 
eraged over time and over 10 different simulations; it is 
also normalized to the external applied strain, in order to 
quantify the discrepancy with respect to a laminar flow. 
The results, as a function of the applied external force 
F , are plotted in Fig. [TT] 

For each pattern we see that, for sufficiently high values 
of Fq, the flow is laminar, that is the ratio < 7 > /j ex t 
tends to unity. For smaller values of the drag forces, the 
behavior of the systems is rather different and strictly de- 
pendent on the equilibrium configuration. In particular 
the cluster phase is easy to shear. In this case the sys- 
tem is never jammed for the values of Fq we have used, 
suggesting that if a yield stress exists, it is very small. 
When the system is already arranged on stripes parallel 
to the walls, the flow is laminar even for the lowest values 
of Fq. In the case of perpendicular stripes the strain rate 
profile is characterized at first by a plateau, which is con- 
nected to the coexistence of the flowing stripes nearby the 
walls and the jammed central stripes. Then at Fq > 0.03 
the stripes breaks and tend to rearrange in the flow di- 
rection. Finally in the bubble phase, the system is at 
first completely jammed as it demonstrated by the zero 
value of the ratio < 7 > / jext- It seems however that in 
case of perpendicular stripes, the partial jamming per- 
sists on a wider range of Fo with respect to the bubble 
phase. Finally we observe that in the cases in which the 
system undergoes a structural change, such as that from 
the perpendicular alignment of stripes to the parallel one, 
or from the melting of the bubble phase to reform paral- 
lel stripes, a peak occurs in the strain-stress curve. The 
height of such a peak in the bubble phase is greater than 
the corresponding one in the stripe phase for an order of 
magnitude. Moreover strong fluctuations appear in the 
stress-strain curve, which might suggest similarities with 
systems in which avalanche-like rearragements occur. 
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FIG. 11: (Color online)The average strain rate, divided by the 
external strain rate ■yext = Vo/(2L x ), as a function of the exter- 
nal force Fo for different pattern morphology at equilibrium. The 
cluster case corresponds to p = 0.15 L x = 76, the parallel stripes 
to p = 0.3 L x = 53, the perpendicular stripes to p = 0.3 L x = 48, 
the bubble case to p = 0.5 L x = 76. 

parallel to the flow direction. This fact is particularly 
evident at high densities (per 2 > 0.3). In particular the 
parallel orientation appears to be preferred even when 
the geometrical constraints impose an equilibrium per- 
pendicular orientation or when the equilibrium pattern 
has a completely different symmetry such as in the bubble 
phase at high density. The rheological behavior discussed 
here has been obtained at T = 0. It is thus important 
to discuss which features are expected to persist when 
thermal effects are taken into account. We can expect 
that if the temperature is below the melting transition, 
the behavior observed at high shear rates should be inde- 
pendent on the temperature. At low shear rates, close to 
the jamming transition, thermal fluctuations could help 
overcome the geometrical contraints leading to a slow 
creep deformation. It would also be interesting to study 
the response of the model under the combined effect of 
thermal fluctuations and oscillating shear stresses. 



The main conclusion of our study is that for sufficiently 
strong shear rate the system dynamically forms stripes 
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